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Abstract 

This report describes the development and validation of LINFLUX-AE, a turbomachinery aeroelastic 
code based on the linearized unsteady 3-D Euler solver, L1NFLUX. A helical fan with flat plate geometry 
is selected as the test case for numerical validation. The steady solution required by L1NFLUX is obtained 
from the nonlinear Euler/Navier Stokes solver TURBO-AE. The report briefly describes the salient features 
of L1NFLUX and the details of the aeroelastic extension. The aeroelastic formulation is based on a modal 
approach. An eigenvalue formulation is used for flutter analysis. The unsteady aerodynamic forces required 
for flutter are obtained by running L1NFLUX for each mode, interblade phase angle and frequency of 
interest. The unsteady aerodynamic forces for forced response analysis are obtained from L1NFLUX for the 
prescribed excitation, interblade phase angle, and frequency. The forced response amplitude is calculated 
from the modal summation of the generalized displacements. The unsteady pressures, work done per cycle, 
eigenvalues and forced response amplitudes obtained from L1NFLUX are compared with those obtained 
from L1NSUB, TURBO-AE, ASTROP2, and ANSYS. 


Introduction 

An overview of the aeroelastic analysis methods for turbomachines (ref. 1 ) shows both time and 
frequency domain methods have been used to obtain unsteady aerodynamic forces and to solve the 
aeroelastic equations. It was noted that time domain methods require large computational times compared to 
frequency domain methods, and should only be used when nonlinearities are expected, and when the need 
justifies the cost. Two approaches were used in obtaining the unsteady aerodynamic forces using frequency 
domain methods. In the first approach (ref. 2) the unsteady aerodynamic equations are linearized about a 
uniform steady flow, thereby neglecting the effects of airfoil shape, incidence and thickness. The unsteady 
aerodynamic models developed in references 3 to 6 based on this approach were used in references 7 and 8 
to study the flutter and forced response analysis of a compressor rotor and propfans. However, methods 
developed by the uniform steady flow approach are restricted to shock-free flows through lightly loaded 
blade rows. 

In the second approach (ref. 9) the unsteady non-linear aerodynamic equations are linearized about the 
non-uniform steady flow. This results in variable coefficient linear unsteady aerodynamic equations, which 
include the effects of steady aerodynamic loading due to airfoil shape, thickness and angle of attack. 
Following the second approach, unsteady linearized Euler aerodynamic models that include the effect of 
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steady aerodynamic loading were developed in references 10 to 14. The aeroelastic code described here is 
based upon the two and three-dimensional linearized Euler codes named respectively L1NFLUX2D and 
L1NFLUX, developed in references 13 and 14 under a NASA contract. These codes were based respectively 
on the non-linear Euler solver developed in references 15 to 17. 

In reference 18, the unsteady aerodynamic calculations from L1NFLUX2D were used with MISER 
(ref. 7), an aeroelastic stability and response code based on a typical section structural model. In 
reference 19, the unsteady aerodynamic calculations from L1NFLUX2D were used with ASTROP2 (ref. 8), 
an aeroelastic stability and response code that uses strip theory to integrate two-dimensional aerodynamic 
forces on to a three dimensional structure. Flutter and forced response calculations were presented for 
cascades in subsonic and transonic flow, with and without mistuning. The aeroelastic calculations with 
MISER and ASTRO P2 approximate the behavior of a three dimensional structure with a two dimensional 
equivalent. However, a more accurate aeroelastic analysis would couple a three dimensional structural 
analysis with a three dimensional aerodynamic analysis. 

The primary objective of the present study is to develop a 3-D-aeroelastic code by coupling the 3-D 
linearized Euler aerodynamic analysis code, L1NFLUX of reference 14, with a three-dimensional structural 
dynamic analysis model. A modal approach is used for flutter and forced response analysis. The resulting 
code is called L1NFLUX-AE to indicate the AeroElastic extension to L1NFLUX code. Brief descriptions of 
the formulation and method of analysis are given in the next section, followed by results and concluding 
remarks. 


Formulation 

The aeroelastic formulation involves solutions from a linearized unsteady aerodynamic model and an 
aeroelastic model coupling the unsteady aerodynamic solution with the structural dynamic solution. The 
aerodynamic model and the aeroelastic formulation are described in this section. 


Aerodynamic model 

The unsteady aerodynamic solution is obtained from the linearized Euler solver, L1NFLUX, developed 
in reference 14. In this code, the linearized unsteady Euler equations are obtained by expanding the 
dependent variables in the unsteady non-linear Euler equations in an asymptotic series of the form 

U= U (x) + u(x(x, t ), t) + higher order terms (1) 

Where x is the mean position, and x is the instantaneous position. Assuming the unsteady excitations 
are harmonic in time, and with the first order variable to be represented as complex valued, the above 
equation can be written as 


U = U(x) + Re{«(x)exp(/co?)} (2) 

Here, the term U (x) is of order one and the second term is of the order £. 

Substituting the expansion of equation (2) in the non-linear unsteady Euler equations, and equating 
terms of like power in £, and neglecting terms of second and higher order in £, nonlinear steady equations 
and linear variable coefficient unsteady equations are obtained. 

For harmonic blade motions with constant phase angle between adjacent blades (interblade phase 
angle), the values of interblade phase angle (a) that can occur are given as (ref. 21). 
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<3,. = 2 w/ N ; r = 0,1,2, 


,N- 1 


( 3 ) 


where X is the number of blades in the blade row. In a time domain approach with periodic boundary 
conditions, the number of blade passages required to be modeled depends on the interblade phase angle, and 
small phase angles may require large number of blade passages to calculate the unsteady aerodynamic 
forces. However, for the linearized approach with assumed harmonic variation in time, the periodic 
conditions are applied on a single extended blade passage region i.e., a region of angular pitch, 

Q = 2n/N (4) 

In solving the linearized unsteady equations, the dependent variables are regarded as pseudo-time 
dependent. This allows solutions to be determined using conventional time-marching algorithms to 
converge the steady and the complex amplitudes of the unsteady conservation variables to their steady state 
values. For more details, see reference 14. 

The steady aerodynamic loading required for the linearized solver, L1NFLUX, is obtained from the 3-D 
Euler/Navier-Stokes solver TURBO-AE, (ref. 20). The code is based on a finite volume scheme. Flux 
vector splitting is used to evaluate the flux Jacobians on the left-hand side of the governing equations and 
Roe's flux difference splitting is used to form a higher-order Total Variation Diminishing (TVD) scheme to 
evaluate the fluxes on the right hand side. Newton sub-iterations are used at each time step to maintain 
higher accuracy. A Baldwin-Lomax algebraic turbulence model is used in the code. Additional details 
regarding the code are available in references 16, 17, and 20. Only the Euler option of the TURBO-AE code 
is used in the present study. 


Aeroelastic model 

The aeroelastic equations of motion for each blade of the blade row of a turbomachine can be written as 

[m\x}+ [k]{x}= {P(X, t)}+ {F(0} (5) 

where [M\ is the mass matrix, [A] is the stiffness matrix that includes the effects of rotation (centrifugal 
stiffening and softening), X] is the vector of blade deflections at the finite element grid points, {P(X,t)} is 
aerodynamic force vector due to blade vibration, and F(t ) is the aerodynamic force vector due to unsteady 
gust or excitation independent of blade vibration. The motion dependent forces cause flutter, and motion 
independent forces cause forced response (forced vibration). 

A modal formulation is used to proceed further. First, a free vibration solution of the blades giving 
generalized masses, natural frequencies, and mode shapes is obtained. The general vibratory motion can 


then be expressed as a superposition of the contributions of the normal modes, [())] = [c|)| , 4> 2 , Anm ] > as 

i=NM 

( 6 ) 
/■= i 


where {q} is the generalized displacement, and NM is the number of normal modes to be used in the 
analysis. Substituting equation (6) in equation (5), and post multiplying the result by [O] , the equations of 
motion for the s th blade can be written as 

[M s ]fe, }+ Us fo s }= [A s hs }+ {AD S } (7.1) 
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where 


\M s ]=[*Y\MMl 

(7.2) 

UJ = [#TO)j}; 

{^}=[<|>F{f(0} 

[M s ] and [AT,] are generalized mass and stiffness matrices, which are diagonal, {q s }is the generalized 
displacement vector, [A} is the blade vibration dependent generalized aerodynamic load matrix, and {AD S } 
is the blade vibration independent generalized aerodynamic load vector. The matrices [MJ, [ AT V ] , and [A s ] 
are of NM x NM size, {q s } and {AD S } are of NM x 1 size. 

For a tuned cascade (or rotor), in which all the blades are identical, the aeroelastic modes consist of 
individual blades vibrating with equal amplitudes at a fixed interblade phase angle between adjacent blades. 
Hence, the motion of the s th blade in r th interblade phase angle mode can be written as 

k } = [q os v “ = few y° ; - v“ (8) 


where c r is given by equation (3). 

Thus the equations of motion for the s th blade, equation (7.1) becomes 

- CO 2 [M s \q ar }e l{at+a A) + [K s \q ar = [ Ay (®+<v) + y(m + a r s) (9) 


Since the blades are identical, the same equation is obtained for each blade. Thus, no additional information 
can be obtained by assembling the equations for all the blades on the disk. Instead, equation (9) can be 
solved for N different values of the interblade phase angles given by equation (3). It should be noted that for 
a mistuned rotor, where the blades have different structural and/or aerodynamic properties, equation (7) 
have to be assembled for all the blades in the rotor and solved, see reference 7 for more details for mistuned 
formulation. 

Calculation of elements of[M s ], [A/], and [<()]. — The elements of [M s ], [A/ v ] and [4>] can be obtained by 
the free vibration analysis using a commercial code like ANSYS/ NASTRAN or by writing an analysis 
program. The matrices [M s ] and [A/] are diagonal. The elements of [A7 V ] and [A/] are related as 

Ku=M u cok (10) 

where co„ is the natural frequency of the i th mode. It is to be noted that usually the mode shapes are mass 
normalized giving M„ =1.0. 

Calculation of elements of [A,] and {AD ,}. — The LINFLUX code is used to obtain the unsteady 
pressures required for calculating the elements of [ A ,. ] and {AD,} matrices. For flutter, the analysis is carried 
out for a given frequency and interblade phase angle, V, and for a given mode of vibration. The output is 
the real and imaginary components of unsteady pressures. This is repeated for N interblade phase angles 
given by equation (3), for NM modes of vibration, and for a range of frequencies of interest. In the case of 
forced response for an external excitation, the analysis is carried out for the unsteady excitation at the 
required interblade phase angle and frequency of interest. 

Calculation of elements of [A,]: The elements of [A,] are due to the vibration of the blade in a given 
mode. Once the unsteady aerodynamic pressures are calculated using LINFLUX for a given frequency and 
interblade phase angle, the elements of [ A r ] are given as 
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A nj =\(>r dA Prj 


( 11 . 1 ) 


where p rj is the unsteady pressure due to blade vibration in j th mode for the r th interblade phase angle 
obtained from LINFLUX, 5, is the i th modal deflection, and dA is the elemental area. The elements of 
equation (1 1 . 1 ) are given by 


n=nnodes 

A rij ~ * Prjn * dA n (11.2) 

n= 1 

where ‘ nnodes ’ is the number of nodes on the aerogrid. Equation (1 1.2) can be expanded as 

n=nnodes 

A rij ~ 'Y, Prjn v^axin * dAx n +& a yi n *dAy n + 5 az j n * dAz n ) (11.3) 

n- 1 

where S UT „ 8 flV , and 8 az , are the modal values at the node n for the i th mode, and dAx, dAy and dAz are the area 
contributions to the node n. 

It should be noted that the pressure in LINFLUX is normalized by (p * aj), where p is the reference air 
density and a„ is the reference speed of sound. In calculating area and modal values, attention should be 
given to non-dimensionlization of length scaling. If L is the reference length used for scaling the geometry 
and the modal values, the final A rij have to be multiplied by L**4 (L**2 for modal values and L** 2 for 
area). It should also be noted that the modal values (mode shape) should be divided by the reference length 
when input to LINFLUX. 

Calculation of elements of {AD r } : The elements of {AD,} are due to the external unsteady disturbances 
coming onto the blade and may be due to acoustic, entropic, vortical or wake interactions. These 
disturbances are independent of the blade vibratory motion, and therefore contribute only to the forced 
response aspect of the aeroelastic problem. The LINFLUX solver calculates the unsteady pressures due to 
these excitations for a given interblade phase angle, r. Let L p) be the unsteady pressure due to these 
excitations on the blade surface. Then the elements of {AD,) are given by 

AD ri =\l t .dA Pr (12.1) 

where 8, is the i th mode shape. The elements of equation (12. 1) are given by 

n=nnodes 

AD ri = Y ^ in * dd n * Pm ( 1 2 - 2 ) 

n - 1 


where ‘ nnodes ’ is the number of nodes on the aerogrid. Equation (12.2) can be written as 


n=nnodes 


AAA ri ~ Y^Pm^axin *dAx n + § a yi n *dAy n +§ az i n * dAz n ) 


(12.3) 


n = 1 

where 8 UT „ 8 flV ,-, and 8 az , are the modal values at the node n for the i th mode, and dAx, dAy and dAz are the 
area contributions to the node n. 
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As mentioned before, in calculating area and modal values, attention should be given to non- 
dimensionlization of length scaling. If L is the reference length used for scaling the geometry and the modal 
values, the final AD ri have to be multiplied by L**3 (L for modal values and L**2 for area). 

It should be noted that the mode shape values are interpolated onto the aerogrid before running 
L1NLFUX. At each CFD grid point on the blade surface, the distance to the nearest three finite-element grid 
points is calculated. Then, the modal deflections, [<])], at these three nearest neighbors are used in a bi-linear 
interpolation scheme to calculate the interpolated value of the modal deflection at that CFD grid point. This 
is done considering the blade’s undeflected position as the reference position. The interpolated modal 
deflections 5 are stored and are used by L1NFLUX. 

Stability calculations. 

Work per cycle approach : To determine aeroelastic (flutter) stability using work per cycle approach, the 
blade is forced to undergo a harmonic vibration in a normal mode with a specified frequency. The vibration 
frequency is typically the natural frequency for the mode of interest, but some other frequency can also be 
used. The aerodynamic forces acting on the vibrating blade and the work done by these forces on the 
vibrating blade during a cycle of vibration are calculated as follows. 

The blade motion is specified to be harmonic as 

q(t) = c/o sin(cot) (13) 

where q 0 is the amplitude of motion and CO is the vibration frequency. 

The work-per-cycle, W, done on the blade is calculated as: 

1T = |j" pdA ,{c)X /c)t]dt (14) 

5 

or, 

W = || pdA.5q 0 (i)cos((i)t)clt (15) 

5 

where, p =p(x,y,z,t) is the unsteady pressure on the blade surface due to blade vibration, A is the blade 
surface area vector pointing into the blade surface, J is the integral over the blade surface, j> is the 

5 

integral over one cycle of blade vibration. The equation for work per cycle using the linearized unsteady 
aerodynamic equations is given in reference 14. 

The work-per-cycle is an indicator of aeroelastic stability. The blade is dynamically unstable if the work 
done on the blade during a cycle of blade vibration is positive. In other words, when W< 0 energy is 
dissipated by the blade and the system is stable, but when W> 0, energy is gained by the blade and it is 
unstable. Note that coupled mode flutter cannot be modeled with this approach, since the flutter mode and 
frequency have to be known a priori for calculating the work. 

Eigenvalue approach : An eigenvalue approach allows one to investigate coupled mode flutter. The 
flutter frequencies and flutter modes are calculated rather than assumed as in aerodynamic work approach. 
For a stability calculation (flutter), the motion-independent forces {AD,) are set to zero in equation (9). 
Dropping the subscript s, since each blade is identical, and canceling out the exponential terms, equation (9) 
can be written as 


- co 2 [M]{7}+ [K]{Y}~ [A r ]{Y}= {0} 


(16) 
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where {7} is same as {q ar \. 

Dividing equation (16) with square of an assumed frequency, co 0 (for which the unsteady aerodynamic 
forces are calculated) 


-(co/co 0 ) 2 [m]{7}+M-[4.]{7}/co 0 2 ={0} (17) 

Rearranging, the equation (17) can be written in the standard eigenvalue problem as: 

[fo]-y[0]]fr-}={o} 08) 

where 

[Prhm-UAl^o 2 (19) 

\Q r ]=[M] (20) 

and 

y = (co/co 0 ) 2 (21) 

The solution of the above eigenvalue problem results in NM complex eigenvalues of the form 

i —=iJy= (I + iv (22) 

C0 o 

The real part of the eigenvalue (jl) represents the damping ratio, and the imaginary part (v) represents the 
damped frequency; flutter occurs if (I > 0 for any of the eigenvalues. 

For the tuned cascade, the stability of each phase angle mode is examined separately. The interblade 
phase angle is fixed at one of the values given by equation (3), and the NM x NM eigenvalue problem is 
solved. The value of interblade phase angle is then changed, and the procedure is repeated for each of the N 
permissible values. The critical phase angle is identified as the one, which results in the lowest flutter speed. 

Forced response calculations'. The equations of motion for the response of the blade are obtained from 
equation (9), by setting the motion dependent unsteady aerodynamic forces, [A,.], to zero, as 

~m 2 [M s lq ar }e i{M+a '^ + = {AD r \e i ^ m+ ° rS ^ (23) 

Since the blades are identical, the same equation is obtained for each blade. Thus, no additional 
information can be obtained by assembling the equations for all the blades on the disk. Instead, 
equation (23) can be solved once for the given excitation. 

Dropping the subscript s, since each blade is identical, and canceling out the exponential terms, 
equation (23) can be written as 


-CO 2 [MtY}+[K\Y}={AD r } (24) 

where {7} is same as {q ar } and CO is the excitation frequency. 

The aeroelastic response of the blades is calculated from equation (24) as 
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(25) 


{7 r }=[M-co 2 [M]]" 1 {^} 

The amplitude of the blade, {X}, is obtained by summing contributions of {7} from all the modes. 

i=NM 

{X}=[#{7}= (26) 

i = 1 


Procedure for Aeroelastic Analysis 

The procedures to do flutter analysis and forced response analysis are shown in figures 1.1 and 1.2 
respectively. The procedures consist of running six programs: (1) TURBO-AE, (2) INTERFACE, (3) a 
structural dynamic analysis solver and the corresponding processing code, RDVIB, (4) a PRE-processor, (5) 
LINFLUX, and (6) a POST-processor. The function of each program is explained in table 1. The response 
analysis differs from flutter analysis in that for flutter analysis the LINFLUX code has to be run for a 
selected number of modes, frequencies, and phase angles, whereas for response analysis, the LINFLUX 
code has to be run only once for the given excitation. 

The description of the steps to obtain an aeroelastic solution is as follows: In step 1, a steady 
aerodynamic solution is obtained from TURBO-AE for the blade. The steady aerodynamic solution is 
written as a database. Step 2 consists of running the INTERFACE program to rewrite the steady data base 
to the format required by LINFLUX. In step 3 a vibration analysis is carried out for the blade to obtain 
structural grid coordinates, generalized masses, natural frequencies, and mode shapes. The code RDVIB is 
used to convert these outputs to a format required by the PRE-processor. In step 4 the PRE-processor is run 
to interpolate modal values at the structural grid nodes onto the aero grid nodes, for the required number of 
modes. In step 5, the LINFLUX code is run for the mode of interest and frequency. At this stage LINFLUX 
also calculates work per cycle for the mode of interest. In step 6 the POST processor is run, to calculate the 
generalized forces, to calculate flutter using the eigenvalue approach, and to calculate the forced response 
amplitudes. 


Results 

In this report, the verification of unsteady pressures, work done per cycle, eigenvalues, and forced 
response amplitudes is given for a helical fan with flat plate geometry. The unsteady pressures calculated 
using the LINFLUX code are compared with those obtained from the LINSUB code (ref. 2). The LINSUB 
code is a two-dimensional unsteady cascade aerodynamic code based on linear unsteady aerodynamic 
theory of Smith (ref. 3). The LINFLUX calculated work done per cycle is compared with those obtained 
from LINSUB and TURBO-AE codes. The TURBO-AE code is a three-dimensional Euler aeroelastic code 
developed in reference 20. The eigenvalue calculations are compared with those obtained from the 
ASTROP2 code of reference 8. The ASTROP2 code uses strip theory to integrate two-dimensional 
aerodynamic forces on to a three dimensional structure. It uses an eigenvalue approach to calculate flutter 
stability. The LINFLUX calculated displacement amplitudes are compared with those obtained from the 
AN SYS finite element program. 
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Helical Fan Geometry and the Aero Grid 


The helical fan configuration consists of a rotor with twisted flat plate blades enclosed in a cylindrical 
duct with no tip gap. This configuration was developed by researchers to provide a relatively simple test 
case for comparison with two-dimensional analyses. Note that there is no experimental data available for 
this configuration. 

The parameters for this configuration are such that the mid-span location corresponds to a flat plate 
cascade with a stagger angle of 45°, unit gap-to-chord ratio, operating in a uniform mean flow at a Mach 
number of 0.7 parallel to the blades. The rotor has 24 blades with a hub/tip ratio of 0.8. The radius at the 
hub is 8.619 cm (3.395 in.) and the radius at the tip is 10.775 cm (4.244 in.). The inlet flow (axial) Mach 
number is 0.495, and the rotation speed of the fan is 16,962.4 rpm, giving a relative Mach number of 
approximately 0.7 at the mid span section. 

A 141 by 11 by 41 grid is used for the calculations. On each blade surface, 81 points are located in the 
chordwise direction, 1 1 points in the spanwise direction, and 41 points in blade to blade direction. The inlet 
and exit boundaries are located at an axial distance of approximately 0.7 chord lengths from the blade 
leading and trailing edges. 


Steady Aerodynamic Solution 

A steady flow solution is obtained first with the TURBO-AE solver. The input to TURBO-AE consists 
of free stream Mach number, rotational speed, time step/CFL number, and back pressure ratio. The code is 
run until the residuals decay to a very low order. For the present case, a time step of 0.04488 is used giving 
a max CFL number of 60.0. About 1000 steps were needed to obtain the desired convergence. Reference 14 
suggests that time discretization errors may result in spurious entropic / vortical boundary layer on the 
surface, which may not allow for the L1NFLUX code be run for unsteady solutions. It is recommended that 
for blade motion and vortical excitations the total pressure loss on the blade surface should be less than 1 0 
percent. 

At the end of the execution of TURBO-AE, the steady aerodynamic solution and input data are written 
to a file. The program INTERFACE is run next to convert the steady aerodynamic solution data obtained 
from TURBO-AE to the form required by L1NFLUX. 


Validation of LINFLUX for Unsteady Aerodynamic Solutions 

Unsteady solutions for prescribed blade vibrations: internally calculated modal values. — To validate 
the LINFLUX code for unsteady solutions, unsteady pressures are calculated for harmonic blade vibration 
in a prescribed mode. The prescribed mode shapes are rigid plunging motion perpendicular to the chord and 
rigid pitching motion about the leading edge. The amplitude of vibration does not vary along the span. This 
choice of mode shapes is meant to reduce the three dimensionality of the unsteady flow field for ease of 
comparison with two-dimensional analyses. 

LINFLUX can be run by an input choice, iflut = 1 , by having these prescribed mode shapes internally 
calculated. However, the purpose of this effort is to run for arbitrary mode shapes that are calculated outside 
LINFLUX and to see that they are correctly read by LINFLUX. To check this, the mode shapes are 
obtained by running LINFLUX with iflut = 1 and are written on unit 3. PRE is then run to interpolate the 
mode shapes onto the aerogrid. But for this case this will not change the modal values at the aerogrid points 
since the modal values are obtained from aerogrid to start with. So actually step 4 is not required for this 
exercise. After running PRE the interpolated mode shapes are written to a file. Now LINFLUX is run with 
iflut = 2, which reads the externally supplied mode shape values at aero-grid points. LINFLUX is run for 
unsteady pressures and compared with those obtained with iflut = 1 and published results. 
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Figures 2 and 3 show the unsteady pressures at mid span obtained for a reduced frequency, CO, of 1.0 
(co = co 0 LIU= (f) 0 L/ a oa M r = co/ M r , where cb = (j) 0 L! , and L is the reference length). It should be 
noted that cb is the input to L1NFLUX. For the present calculations the value of speed of sound, <x„, is 
assumed as 13707 in./sec; reference length is the chord which is 1.0 in., giving a value of 1527 cycles/sec 
for the assumed frequency co 0 . Figure 2 shows the variation of unsteady pressure with chord for plunging 
motion and Fig 3 for pitching motion about leading edge. Fig 2. 1 is for zero phase angle, and Fig. 2.2 is for 
180° phase angle. Figures 3.1 and 3.2 show the unsteady pressure distributions for pitching motion for zero 
and 90° phase angle respectively. The pressures are compared with those obtained from L1NSUB code 
based on Smith theory (ref. 3). To compare with L1NSUB results, the unsteady pressure output from 
L1NFLUX is divided by square of the Mach number times the amplitude of oscillation. The pressure 
distribution predictions from LINFLUX show excellent agreement with those from L1NSUB. The work 
done per cycle calculations, not shown here, also showed excellent agreement with L1NSUB predictions. 

Unsteady solutions for prescribed blade vibrations: externally supplied mode shapes. — Further 
verification of iflut = 2 option carried out by creating a modal file outside the LINFLUX code. A program is 
written to simulate rigid plunging motion and rigid pitching motion about mid chord for the helical fan. The 
amplitude of pitching is 1 percent chord and pitching is 2° (0.03491 radians). The modal shape values thus 
generated are written to a file. LINFLUX is run for these externally supplied mode shapes using iflut = 2 
option. Figures 4 and 5 show the unsteady pressure difference at mid span. At this section co = 1.0, (O = 0.7 
and the Mach number M = 0.7. Figures 4.1 and 4.2 show plunging motion results for ct = 0 and 180°, and 
figures 5.1 and 5.2 show results for pitching motion about mid chord for ct = 0 and 180°. They are 
compared with those obtained from LIN SUB. The LINFLUX calculations show excellent agreement with 
L1NSUB calculations for both plunging and pitching motion, since the simulated motion is close to a two 
dimensional rigid plunging and pitching motion. The work done per cycle obtained from LINFLUX, 
L1NSUB and TURBO-AE is compared in figures 6.1 and 6.2 for pitching and plunging motions 
respectively. The LINFLUX calculations show excellent agreement with L1NSUB predictions. However, 
TURBO-AE predictions show a slight difference from L1NSUB calculations. This may be because in 
TURBO-AE reflecting boundary conditions are used in obtaining the solutions. Table 2 shows the 
numerical values of the work done per cycle (wpc) for both the modes from LINFLUX, TURBO-AE and 
L1NSUB. It should be noted that the LINSUB wpc values are multiplied with (amplitude of oscillation**2) 
(M**2)*(span) with span equal to 0.849 for the present case. The excellent agreement in unsteady pressures 
validates the unsteady pressure calculations with LINFLUX for externally supplied modes. 

Unsteady solutions for external excitations. — To validate LINFLUX for forced response excitation, 
analysis is carried out for three types of excitations: a downstream acoustic excitation, upstream acoustic 
excitation, and a vortical gust. The unsteady pressure difference is compared with those obtained by 
LINSUB in figures 7 to 9. It should be noted that structural mode shapes are not required for running 
LINFLUX for external excitations. 

Figures 7 and 8 show the unsteady pressure difference calculated for to = 3.333 for upstream and 
downstream acoustic waves for a = 90° respectively. The results for the upstream acoustic wave show 
excellent agreement with those obtained from LINSUB. The results for downstream acoustic wave show 
fair agreement with those obtained from LINSUB. In the case of upstream acoustic wave the wave velocity 
is same as the free stream velocity, and the grid resolution used in front of the cascade is sufficient to 
resolve the wave velocity. However, the downstream acoustic wave is traveling opposite to the free stream, 
and the grid aft of the cascade may require higher resolution to match the LINSUB predictions. 

Figure 9.1 shows the unsteady pressures for a vortical gust at -90° phase angle and CO = 1.287, and 
figure 9.2 shows the unsteady pressures for vortical gust at -180° phase angle and co = 2.574. The unsteady 
pressures are compared with those obtained from LINSUB. The figures show that the trends are correctly 
predicted but the magnitudes are not predicted well. This may indicate that for the vortical gust excitations a 
finer grid is required to simulate the wave lengths correctly. 
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Structural model, mode shapes, and frequencies 


The purpose of this development is to use L1NFLUX for aeroelastic analysis of real blades with real 
mode shapes with a mode shape file that is externally supplied to LINFLUX. For this purpose a realistic 
structural dynamic analysis for structural mode shapes, generalized masses, and frequencies is required. 

This will help in validating the modules PRE and POST, and streamlining the procedure for flutter and 
forced response analysis of real blade geometries. For the helical fan geometry considered here there was no 
structural model available. Therefore, from the available aerodynamic grid, a finite element grid was made. 
The finite element grid is made exactly to match with the aerogrid by using 8 1 nodes along the chord and 
1 1 nodes along the span giving 891 nodes on the surface. The blade is assumed to have thickness of 0.02 
inches. The material of the blade is assumed to be titanium with a material density of 0.160 lb/in**3, 
Young’s modulus is 16xlOE6 lb/in**2 and the Poisson’s ratio is 0.3. The blade is rotating at 16,962 rpm 
about the engine axis. 

The structural dynamic analysis was performed with the AN SYS computer code. SFIELL 63 elements 
were used in the analysis. Dynamic analysis included the effects of spinning, and stress softening. Figures 
10.1 and 10.2 show the blade modes shapes for first and second modes. The first mode is a bending mode at 
822 FIz and the second mode is a torsion mode at 1882 FIz. The frequency values are tabulated in table 3 for 
both rotating and non-rotating blades. For the present aeroelastic analysis the first two modes are used. 


Validation of PRE 

The purpose of the module PRE is to rotate the structural gird to aero grid coordinate system, and to 
interpolate modal values from structural grid onto aerogrid. These interpolated modal values are required to 
calculate the unsteady aerodynamic forces using LINFLUX and to calculate the generalized forces in 
POST. Any non-dimensionalization of the geometry is also to be accounted for at this stage. After running 
PRE the interpolated mode shapes are written to a file. The interpolated modal values are visually checked, 
and also by plotting the interpolated mode shapes. Table 4 shows the modal values on the structural grid 
and the interpolated values on the aerogrid at selected nodes. They show excellent agreement as expected 
since the structural and aerogrid are identical for the present example. 


Validation of POST 

In POST, the unsteady pressures calculated by LINFLUX, and the interpolated mode shapes from PRE 
are used to calculate generalized forces and displacements, eigenvalues, and response amplitudes. The 
following describes how they were systematically checked. 

Generalized forces. — The generalized force calculation is checked by assuming uniform pressure 
on the aerogrid and assuming a modal value of 1.414 perpendicular to chord at all points. The blade face 
area ofthe helical fan is [length * chord = (4.244 - 3.395)* 1 .0] = 0.849 in**2 (calculated area = 0.8491). 

As the helical fan has 45° twist, the area "Ldax = 0.849*0.707 = 0.62923 (0.5998) and 
'Lday = 0.849*0.707 = 0.62923 (0.5913) and 'Zdaz = 0.0 (0.09938), a quantity an order less than 'Ldax and 
'Lday. The displacement values in x,y, and z directions are 1.00, 1.00 and 0.0. Substituting these values in 
equation (1 1.3) one obtains a value of (0.62923*1.0 + 0.62923*1.0 + 0.09938* 0.0) = 1.26 (1.19) for 
generalized force. The computer calculated generalized force values show a 6.0 percent difference from 
hand calculated values. This is expected since the hand calculations assumed a constant twist of 45° along 
the entire span. 

The generalized forces and the displacements calculations are also checked by assuming uniform 
pressure on the aerogrid but using the AN SYS calculated modal values for two modes of the non-rotating 
blade. Both ANSYS and POST use harmonic response analysis method to determine the steady-state 
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response of a linear structure to loads that vary sinusoidally with time. The mode supeiposition method is 
used to calculate the structure’s response. Table 5 shows th ex,y, and z components of the displacements at 
the tip leading edge (node number 8 1 1 ), mid chord (851) and at the trailing edge (891) obtained from both 
POST and ANSYS. The excitation frequency is 400 Hz. They show excellent agreement. 

Flutter . — Flutter analysis with L1NFLUX is carried out using the first two modes given by ANSYS for 
a frequency of 1527Hz ( cd = 0.7 and CO = 1.0). The calculated unsteady pressures are used by POST to 
calculate generalized forces and eigenvalues. To validate the eigenvalue calculations, the present results are 
compared with those obtained from the ASTROP2 code. As mentioned before, the ASTROP2 code uses 
strip theory to integrate two-dimensional aerodynamic forces on to a three dimensional structure. The 
aerodynamic forces are calculated about leading edge, using theory of reference 4. It also uses eigenvalue 
approach to calculate stability. 

Unsteady pressure comparison for the calculated structural mode shapes : Figures 1 1 and 12 show the 
unsteady pressures calculated for model and mode 2 respectively, for a = 0 and 180° at cd = 0.7. They are 
compared with TURBO-AE. They agree qualitatively, but show some quantitative differences. In the 
previous sections it was shown that calculations from L1NFLUX compared well with LIN SUB predictions 
for assumed modes. Since the airfoil geometry is same except for the mode shapes the L1NFLUX 
calculations are assumed to be accurate. The quantitative differences between TURBO-AE and L1NFLUX- 
AE results indicate that a more accurate TURBO-AE solution is required by optimizing the time step, 
amplitude of oscillation, and number of cycles used to obtain the solution. For the present report time step, 
amplitude of oscillation and number of cycles are fixed, which may affect the solution for phase angles, and 
frequencies. In addition, reflecting boundary conditions are used in the TURBO-AE runs, whereas 
L1NFLUX uses non-reflecting boundary conditions. It should be noted that getting an optimum TURBO- 
AE solution is not the objective of this report. 

Work per cycle (wpc): As mentioned before, work per cycle can be calculated with L1NFLUX for each 
mode at a given frequency. This will give information on the stability of each mode (but not for a coupled 
mode for which a coupled mode analysis has to be carried out as given in the next section). The calculated 
work per cycle from L1NFLUX and TURBO-AE are listed in table 6, and shown in figure 13.1 for the first 
ANSYS mode and in figure 13.2 for the second ANSYS mode. All the calculations for 
cd = 0.7. It can be seen that they compare very well. The discrepancies seen are due to the non-reflecting 
boundary conditions used in L1NFLUX and not in TURBO-AE runs. 

Stability root locus: The POST program in L1NFLUX-AE code is used now to calculate the stability 
eigenvalues using the ANSYS calculated mode shapes. The calculations are made for 8 of the 24 possible 
phase angles at cd = 0.7. ASTROP2 calculated cascade parameters and equivalent plunging and pitching 
values at selected strips for this blade are given in appendix A. Table 7 shows the numerical values of 
damping and frequency obtained from L1NFLUX and ASTROP2. Figures 14 and 15 show the plot of the 
frequency versus damping ratio for the two modes. There is only slight difference in the least stable 
damping ratio. This is expected since the L1NFLUX unsteady pressure calculations include the effect of 
chordwise flexibility, and three dimensional flow effects. As expected, the eigenvalues indicate that the 
blades are stable, which was also indicted by the work per cycle approach in the TURBO-AE and 
L1NFLUX programs. These results validate the routines in the post-processor for flutter prediction. 

Forced response amplitudes . — The response amplitudes are checked by using the L1NFLUX calculated 
unsteady pressures for an upstream acoustic wave and using the ANSYS calculated modal values for two 
modes for the non-rotating blades. The frequency of this excitation is 400 Hz. The unsteady pressures from 
the aerogrid are interpolated on to the structural grid and an ANSYS analysis was carried out. Table 8 
shows th ex,y, andz components at nodes 811, 851 and 891 obtained from both POST and ANSYS. They 
show good agreement for both real and imaginary values. 

These results validate the POST module for generalized forces and response amplitude. 
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Concluding Remarks 


The formulation and validation of the turbomachinery aeroelastic code L1NFLUX-AE are presented in 
this report. The code requires execution of six programs: (1) steady aerodynamic solution, (2) interface, (3) 
read structural vibration data, (4) interpolate mode shapes on to aerogrid, (5) a linearized unsteady 
aerodynamic solution, (6) and flutter and response calculations. 

The various modules and the interfaces between the programs have been verified for a helical fan with 
flat plate airfoils. The unsteady pressures calculated from L1NFLUX for both blade vibrations, and external 
disturbances agreed well with those obtained from L1NSUB and TURBO-AE. The flutter eigenvalues from 
L1NFLUX-AE are close to those obtained from ASTROP2. The response amplitudes agreed with those 
obtained by ANSYS harmonic analysis for two modes. 

The helical fan example validated the procedure, and all the programs required in the L1NFLUX-AE 
aeroelastic package for flutter and forced response. Application of the L1NFLUX-AE code for real blade 
geometries is intended in future. 
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Table 1. — LINFLUX-AE programs and their function. 


Number 

Program Name 

Function 

I 

TURBO-AE 

do steady aerodynamic analysis 

II 

INTERFACE 

convert TURBO-AE output to data format required by LINFLUX, write 
on UNIT 51 

III 

RDVIB 

(1) do free vibration analysis. At present done outside this program 
through ANSYS, NASTRAN etc. 

(2) read vibration analysis output file and rotate the grid and mode shapes 
to aero-grid coordinate system, write on UNIT 3 

IV 

PRE 

interpolate (calculate) mode shape values at each aero-grid point on both 
the airfoil surfaces for all modes (UNIT 27). 

V 

LINFLUX 

calculate unsteady pressures for the given mode, frequency and interblade 
phase angle; repeat for all modes, frequencies and interblade phase angles 
(UNIT 94). 

VI 

POST 

calculate generalized forces; flutter eigenvalues and response. 
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Table 2. — Comparison of work per cycle for externally supplied inodes 


(co=1.0, (0 = 0.7, Mat mid section = 0.7). 


Assumed Modes 


LINFLUX 

TURBO- AE 

LINSUB 

Rigid plunging 




0 

-0.16990-03 

-0.17161-03 

-0.17296-03 

45 

-0.22242-03 

-0.20180-03 

-0.22285-03 

90 

-0.26345-03 

-0.22098-03 

-0.25855-03 

180 

-0.46051-03 

-0.41387-03 

-0.46618-03 

225 

-0.38399-03 

-0.34110-03 

-0.38879-03 

270 

-0.28268-03 

-0.22165-03 

-0.28638-03 

315 

-0.22321-03 

-0.21899-03 

-0.22625-03 

Rigid pitching about 
midchord 




0 

-0.54063-03 

-0.53604-03 

-0.57368-03 

45 

-0.43992-03 

-0.29169-03 

-0.46181-03 

90 

-0.54835-03 

-0.52927-03 

-0.55458-03 

180 

-0.13668-02 

-0.15622-02 

-0.13974-02 

225 

-0.14754-02 

-0.14969-02 

-0.15145-02 

270 

-0.14312-02 

-0.11703-02 

-0.14766-02 

315 

-0.12991-02 

-0.10945-02 

-0.13531-02 


LINSUB output is multiplied by amplitude**2*M**2*span of the blade (=0.849). 


Table 3, — Natural frequencies for the helical fan. 


Mode 

F requency, Hz 

n = o.o 

F requency, Hz 
Q = 16962 RPM 

1 

364.0 

822.0 

2 

1745.0 

1882.0 

3 

2220.0 

2517.0 

4 

2538.0 

2992.0 

5 

3667.0 

4046.0 

6 

5075.0 

5165.0 

7 

5472.0 

5803.0 

8 

6507.0 

7530.0 

9 

7499.0 

8176.0 

10 

8321.0 

8877.0 
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Table 4. — Comparison of interpolated modal values. 


Mode 

1 


Lefem 

(811) 

LE Inter 

Mid chord 
(851) 

Inter 

TE 

(891) 

TE inter 


UX 

-757.7949 

-756.6143 

-792.5469 

-790.8140 

-777.0029 

-777.0532 


UY 

-793.9869 

-792.8906 

-837.8209 

-836.0615 

-835.8530 

-835.9058 


UZ 

207.3049 

205.6547 

155.4340 

154.8228 

62.2560 

62.2633 


Mode 

2 


LEfem 

LE Inter 

Mid chord 

Inter 

TE 

TE inter 


UX 

-1483.5699 

-1452.766 

-283.6709 

-277.6899 

1428.260 

1428.260 


UY 

-1332.3199 

-1305.039 

-259.3320 

-253.8341 

1311.449 

1311.443 


UZ 

305.4089 

299.0735 

64.5855 

63.5685 

-136.7429 

-136.7376 


Table 5. — Response amplitudes for uniform pressure at 400 Hz; 


(two modes of the non-rotating blade used) 


Method 

NODE 

UX 

UY 

UZ 


811 

-0.28392 

-0.29496 

0.07645 

POST 

851 

-0.29432 

-0.30894 

0.05721 


891 

-0.28647 

-0.30635 

0.02277 


811 

-0.27556 

-0.28607 

0.07459 

ANSYS 

851 

-0.28451 

-0.29959 

0.05010 


891 

-0.27444 

-0.29379 

0.02181 


Table 6. — Comparison of work per cycle (CO = 1.0, CO =0.7). 



ANSYS Modes 


LINFLUX 

TURBO-AE 

Mode 1: 



0 

-4.9289 

-5.241 

45 

-6.0959 

-5.635 

90 

-7.0073 

-6.743 

135 

-12.601 

-12.78 

180 

-11.986 

-16.00 

225 

-10.343 

-12.49 

270 

-8.0198 

-8.455 

315 

-6.619 

-7.248 

Mode 2: 



0 

-15.053 

-17.85 

45 

-12.226 

-9.283 

90 

-14.383 

-13.93 

135 

-32.38 

-33.28 

180 

-32.66 

-43.89 

225 

-35.55 

-46.80 

270 

-34.821 

-37.14 

315 

-31.902 

-35.66 
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Table 7. — Eigenvalue Comparison co 0 = 1527 Hz, CO = 0.7, a „ = 13707 in /sec,/> 00 = 1 1.881 psi. 


Table 7.1 .— 0 = 0°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.00928 

0.5334 

-0.01426 

0.53448 

2 

-0.05241 

1.17275 

-0.08748 

1.15085 


Table 7.2.— a = 45°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.0116 

0.53408 

-0.01796 

0.53716 

2 

-0.03802 

1.15010 

-0.06256 

1.12055 


Table 7.3.— a = 90°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.01368 

0.53899 

-0.02038 

0.54540 

2 

-0.04464 

1.145 

-0.07099 

1.10848 


Table 7.4.— o= 135°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.02563 

0.54344 

-0.02389 

0.56952 

2 

-0.08755 

1.09407 

-0.21854 

1.0858 


Table 7. 5.— o= 180°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.02415 

0.54218 

-0.03735 

0.56126 

2 

-0.09075 

1.09911 

-0.16117 

1.00596 


Table 7.6.-0 = 225°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.02054 

0.54183 

-0.02917 

0.55654 

2 

-0.10756 

1.11348 

-0.18212 

1.03140 
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Table 7.7 .— 0 = 270°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.01561 

0.53975 

-0.02074 

0.54952 

2 

-0.11459 

1.13596 

-0.19597 

1.07218 


Table 7.8.-0 = 315°. 


Mode 

LINFLUX-AE 

ASTROP2 


Eigenvalue 

Eigenvalue 


real 

imaginary 

real 

imaginary 

1 

-0.01267 

0.53723 

-0.01692 

0.54406 

2 

-0.11217 

1.1596 

-0.2046 

1.11698 


Table 8. — Comparison of response amplitudes for pressure due to upstream acoustic wave 
(00=3.332, o = 90o; two modes of the non-rotating blade used at the excitation frequency = 400 Hz). 


Method 

NODE 

UX 

UY 

UZ 


811 

(-0.34170,-0.19651) 

(-0.35541,-0.20379) 

(0.92197-01,0.52750-01) 

POST 

851 

(-0.35584,-0.20114) 

(-0.37373,-0.21112) 

(0.69180-01,0.39128-01) 


891 

(-0.34941,-0.19239) 

(-0.37340, -0.20624) 

(0.27843-01,0.15262-01) 


811 

(-0.33145,-0.18983) 

(-0.34449,-0.19668) 

(0.89914-01,0.51209-01) 

ANSYS 

851 

(-0.34509,-0.19318) 

(-0.36339, -0.20340) 

(0.60738-01,0.34050-01) 


891 

(-0.33571,-0.18356) 

(-0.35895,-0.19692) 

(0.26731-01,0.14546-01) 
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Figure 1 . 1 . — LINFLUX-AE flow chart for flutter. 
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run TURBO 


: steady aerodynamic solution 



Figure 1.2. — LINFLUX flow chart for forced response. 


NASA/TM— 2004-212978 


20 










Figure 2.1. — Unsteady pressures for plunging motion, 
a = 0°, helical fan, to = 0.7, M = 0.7. 



Figure 2.2. — Unsteady pressures for plunging motion, 
a = 180°, helical fan, to = 0.7, M = 0.7. 
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Figure 3.1. — Unsteady pressure distribution for pitching about 
leading edge, a = 0°, helical fan, to = 0.7 , M = 0.7. 



Figure 3.2. — Unsteady pressure distribution for pitching about 
leading edge, c = 90°, helical fan, to = 0.7, M = 0.7. 
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Figure 4. 1 . — Unsteady pressure comparison for assumed mode, 
plunging, c = 0°, to = 0.7, M = 0.7. 
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Figure 4.2. — Unsteady pressure comparison for assumed 
mode, plunging, o = 180°, to = 0.7, M = 0.7. 
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Figure 5.1. — Unsteady pressure comparison, pitching 
about mid-chord, a = 0°, (0 = 0.7, M = 0.7. 



distance along chord 


Figure 5.2. — Unsteady pressure comparison, pitching about 
mid-chord, a = 180°, di = 0.7, M = 0.7. 
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Figure 6.1. — Work per cycle for assumed plunging mode, helical fan, to = 0.7. 
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Figure 6.2. — Work per cycle for assumed pitching mode, helical fan, di = 0.7. 
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Figure 7. — Unsteady pressure difference, helical fan at r/rd = 0.9, LINFLUX vs LINSUB, 
upstream acoustic wave, kc = 3.332, a = 90°. 



Figure 8. — Unsteady pressure difference, helical fan at r/rd = 0.5, LINFLUX vs LINSUB, 
downstream acoustic wave, (0 = 3.332, <7 = 90°. 
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Figure 9.1. — Unsteady pressures difference, helical fan at r/rd = 0.9, 
LINFLUX vs LINSUB, vortical gust, co= 1.287, a = -90°. 
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Figure 9.2. — Unsteady pressures difference, helical fan at r/rd = 0.9, 
LINFLUX vs LINSUB, vortical gust, oo = 2.574, a = -180°. 
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Figure 10.1. — First mode, 822 Hz. 



Figure 10.2. — Second mode, 1882 Hz. 
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Figure 11.1. — Unsteady pressure distribution for mode 1, helical fan, co = 0.7, a =0°. 
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Figure 1 1.2. — Unsteady pressure distribution for mode 1, helical fan, co = 0.7, a =180°. 
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12. 1. — Unsteady pressure distribution for mode 2, helical fan, to = 0.7, a =0°. 



Figure 12.2. — Unsteady pressure distribution for mode 2, helical fan, to = 0.7, <r =180°. 
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Figure 13.1. — Work per cycle for 1 st ANSYS mode, helical fan, to =0.7. 



Figure 13.2. — Work per cycle for 2 nd ANSYS mode, helical fan, oo = 0.7. 
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Figure 14. — Root locus plot; damping versus frequency for helical fan, co = 0.7, mode 1. 
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Figure 15. — Root locus plot; damping versus frequency for helical fan, co = 0.7, mode 2. 
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Appendix A 

Cascade Geometric and Aerodynamic Parameters 
Calculations for HFAN with ASTROP2 


This appendix presents the calculated values of aerodynamic and geometric parameters obtained from 
the aeroelastic code ASTROP2. These are used in the calculation of flutter eigenvalues in the ASTROP2 
code. The ASTROP2 calculated eigenvalues were compared with the LINFLUX-AE calculations in figures 
15 and 16. The ASTRO P2 code is a quasi-three dimensional aeroelastic code. It combines two-dimensional 
unsteady aerodynamic loads from linear theory of Rao (ref. 4) with the vibration modes and frequencies of 
a three dimensional structure. The aerodynamic loads are combined in a strip wise fashion. An equivalent 
plunging, motion perpendicular to chord, and pitching motion, rotation about leading edge, are calculated at 
a selected number of strips. The unsteady aerodynamic forces are calculated for these plunging and pitching 
motions. 

For the present analysis, nine strips are selected. The location of the strips and their coordinates are 
given in table Al. The length along the blade leading edge is 0.8256 inches. The 75 percent location is at 
3.8990 from the engine axis. 


Table Al. — Strip location and their coordinates. 


Strip number 

Node Number 

X 

Y 

Z 

1 

163 

-0.3654 

3.4449 

0.9256 

2 

244 

-0.3635 

3.5274 

0.9455 

3 

325 

-0.3616 

3.6100 

0.9654 

4 

406 

-0.3597 

3.6925 

0.9851 

5 

487 

-0.3578 

3.7751 

1.0049 

6 

568 

-0.3560 

3.8576 

1.0246 

7 

649 

-0.3542 

3.9402 

1.0442 

8 

730 

-0.3524 

4.0228 

1.0638 

9 

811 

-0.3506 

4.1054 

1.0831 


First step in the ASTROP2 program is to calculate the normal, tangent, and stream vector at the selected 
strips. These values are given in table A2. 


Table A2.1. — Tangent vector. 


Strip number 

Node Number 

X 

Y 

Z 

1 

163 

0.0236 

0.99972 

0.24115 

2 

244 

0.023 

0.99974 

0.24116 

3 

325 

0.023 

0.99974 

0.24032 

4 

406 

0.02303 

0.99973 

0.23943 

5 

487 

0.02238 

0.99975 

0.23931 

6 

568 

0.02178 

0.99976 

0.23717 

7 

649 

0.02178 

0.99976 

0.23729 

8 

730 

0.02178 

0.99976 

0.23729 

9 

811 

0.02178 

0.99976 

0.23181 
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Table A2.2. — Chord line vector. 


Strip number 

Node Number 

X 

Y 

Z 

1 

163 

0.73042 

0.12478 

-0.67151 

2 

244 

0.72263 

0.12416 

-0.67999 

3 

325 

0.71488 

0.12345 

-0.68827 

4 

406 

0.70728 

0.12284 

-0.69618 

5 

487 

0.69969 

0.12212 

-0.70394 

6 

568 

0.69220 

0.12154 

-0.71140 

7 

649 

0.68484 

0.12087 

-0.71860 

8 

730 

0.67876 

0.08017 

-0.72997 

9 

811 

0.67070 

0.00000 

-0.74173 


Table A2.3. — Normal vector. 


Strip number 

Node Number 

X 

Y 

Z 

1 

163 

-0.70141 

0.19199 

-0.72727 

2 

244 

-0.70975 

0.18991 

-0.71959 

3 

325 

-0.71775 

0.18763 

-0.71185 

4 

406 

-0.72540 

0.18538 

-0.70427 

5 

487 

-0.73299 

0.18319 

-0.69678 

6 

568 

-0.74006 

0.17966 

-0.68939 

7 

649 

-0.74711 

0.17816 

-0.68205 

8 

730 

-0.74882 

0.17696 

-0.67685 

9 

811 

-0.74115 

0.17163 

-0.67054 


The calculated stagger, gap to chord ratio, sweep, semi chord , and ship width are given in 
table A3. 


Table A3. — Calculated cascade structural parameters at the strips. 


PT. 

STN 

SWEEP(DEG). 

SEMICHD 

GAP/CHD 

STAGGER 

WIDTH 

1 

3.44490 

1.35205 

0.49609 

0.93961 

41.87744 

0.16995 

2 

3.52740 

1.31790 

0.49873 

0.95678 

42.51256 

0.08489 

3 

3.61000 

1.31790 

0.50143 

0.97365 

43.13840 

0.08498 

4 

3.69250 

1.31968 

0.50412 

0.99040 

43.74040 

0.08484 

5 

3.77510 

1.28207 

0.50687 

1.00681 

44.33916 

0.08496 

6 

3.85760 

1.24804 

0.50969 

1.02290 

44.91832 

0.08484 

7 

3.94020 

1.24804 

0.51252 

1.03881 

45.48364 

0.08491 

8 

4.02280 

1.24804 

0.51513 

1.05438 

46.04495 

0.08491 

9 

4.10540 

1.24804 

0.51997 

1.07228 

47.06836 

0.08484 


The pitching values are calculated by averaging the rotational values along the chord nodes of the strip. 
The effect of sweep is to add contribution due to prime terms in the above table. Since the sweep angle is 
very low the contribution these numbers to aeroelastic stability is negligible. 
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Table A4.1. — Modal values for mode 1. 


PT. 

H-VALUE 

A-VALUE 

HP-VALUE 

AP-VALUE 

HPP- VALUE 

APP-VALUE 

1 

0.701 13E+02 

0.51497E+02 

0.86105E+03 

0.2827 1E+03 

0.47609E+04 

0.59555E+03 

2 

0.15462E+03 

0.67203E+02 

0.1 1638E+04 

0.12313E+03 

0.27826E+04 

0.11671E+04 

3 

0.26084E+03 

0.75581E+02 

0.13728E+04 

0.54733E+02 

0.20162E+04 

0.47436E+02 

4 

0.38294E+03 

0.80152E+02 

0.15292E+04 

0.54810E+02 

0.15581E+04 

0.50465E+02 

5 

0.51702E+03 

0.84879E+02 

0.16489E+04 

0.48828E+02 

0.13008E+04 

0.15708E+03 

6 

0.65975E+03 

0.85496E+02 

0.17402E+04 

0.17189E+02 

0.56030E+03 

0.24749E+03 

7 

0.80822E+03 

0.87237E+02 

0.17557E+04 

0.2365 1E+01 

0.73364E+02 

0.00000E+00 

8 

0.95798E+03 

0.87438E+02 

0.17498E+04 

0.23651E+01 

0.29395E+03 

0.00000E+00 

9 

0.11007E+04 

0.87638E+02 

0.16421E+04 

0.2365 1E+01 

0.18687E+04 

0.00000E+00 


Table A4.2. — Modal values for mode 2. 


PT. 

H-VALUE 

A-VALUE 

HP-VALUE 

AP-VALUE 

HPP-VALUE 

APP-VALUE 

1 

-0.18718E+03 

0.42099E+03 

-0.20246E+04 

0.41826E+04 

-0.52417E+04 

0.83726E+04 

2 

-0.36940E+03 

0.79312E+03 

-0.235 85E+04 

0.47366E+04 

-0.25552E+04 

0.47168E+04 

3 

-0.57641E+03 

0. 12071 E+04 

-0.25654E+04 

0.51218E+04 

-0.20088E+04 

0.40234E+04 

4 

-0.80006E+03 

0.16518E+04 

-0.27344E+04 

0.54204E+04 

-0.19785E+04 

0.2661 1E+04 

5 

-0.10376E+04 

0.21171E+04 

-0.28986E+04 

0.55500E+04 

-0.17507E+04 

0.91943E+03 

6 

-0.12856E+04 

0.259 18E+04 

-0.29430E+04 

0.556 10E+04 

-0.3391 1E+03 

-0.50781E+03 

7 

-0.15370E+04 

0.30539E+04 

-0.29464E+04 

0.53440E+04 

0.41650E+03 

-0.67407E+04 

8 

-0.17783E+04 

0.34785E+04 

-0.27526E+04 

0.47466E+04 

0.34133E+04 

-0.72192E+04 

9 

-0.19931E+04 

0.38538E+04 

-0.23759E+04 

0.41333E+04 

0.55186E+04 

-0.73247E+04 


The blade is rotating counter clock wise at 16,696 lpm about the engine axis. The axial Mach number is 
0.495. The calculated effective Mach numbers, and reduced frequency for an assumed frequency of 
1527.6 Hz is given in table A5. 

REFERENCE FREQUENCY (RAD/SEC & Hz) = 9598.19434 1527.59998 


Table A5. — Calculated ASTROP2 cascade analysis parameters. 


PT 

HMACH 

SWEEP 

EFF. M 

SEMICHD 

RED. FREQ 

GAP/CHD 

STAGGER 

1 

0.67728 

1.35205 

0.67709 

0.49694 

0.51393 

0.93961 

41.87744 

2 

0.68483 

1.31790 

0.68465 

0.49963 

0.51101 

0.95678 

42.51256 

3 

0.69248 

1.31790 

0.69230 

0.50239 

0.50815 

0.97365 

43.13840 

4 

0.70021 

1.31968 

0.70003 

0.50510 

0.50525 

0.99040 

43.74040 

5 

0.70804 

1.28207 

0.70786 

0.50791 

0.50244 

1.00681 

44.33916 

6 

0.71594 

1.24804 

0.71577 

0.51077 

0.49969 

1.02290 

44.91832 

7 

0.72392 

1.24804 

0.72375 

0.51364 

0.49696 

1.03881 

45.48364 

8 

0.73198 

1.24804 

0.73181 

0.51659 

0.49430 

1.05438 

46.04495 

9 

0.74012 

1.24804 

0.73994 

0.51832 

0.49051 

1.07228 

47.06836 
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